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ABSTRACT 

We investigated the dynamical reaction of the central region of galaxies to a falling 
massive black hole by TV-body simulations. As the initial galaxy model, we used an 
isothermal King model and placed a massive black hole at around the half-mass radius 
of the galaxy. We found that the central core of the galaxy is destroyed by the heating 
due to the black hole and that a very weak density cusp (p oc r~ a , with a ~ 0.5) 
is formed around the black hole. This result is consistent with recent observations of 
large elliptical galaxies with Hubble Space Telescope. The velocity of the stars becomes 
tangentially anisotropic in the inner region, while in the outer region the stars have 
radially anisotropic velocity dispersion. The radius of the weak cusp region is larger for 
larger black hole mass. Our result naturally explains the formation of the weak cusp 
found in the previous simulations of galaxy merging, and implies that the weak cusp 
observed in large elliptical galaxies may be formed by the heating process by sinking 
black holes during merging events. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: kinematics and 
dynamics — galaxies: nuclei — galaxies: structure — methods: numerical 
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Introduction 



Recent HST observations ( 


Lauer et al. 1995; 


Byun 


3t al. 1996; Gebhardt et al. 1996; 


Faber et al. 1996|; 


Kormendy et al. 1996 


) have drastically changed our 



understanding of the structure of elliptical galaxies. 
First, no elliptical galaxies have cores with constant 
surface brightness. All of them turned out to have 
cusps with power-law surface brightness profiles. Sec- 
ond, elliptical galaxies are divided into two groups in 
terms of the slope of the cusp: "weak-cusp" galax- 
ies and "strong-cusp" galaxies. The weak-cusp galax- 
ies have the surface brightness profiles expressed as 
I{R) oc Br 1 with <; 7 < 0.3 and deprojected lumi- 
nosity density profiles expressed as p(r) oc r~ a with 



0.3 <J a 1.1 (Kormendy et al. 1996). On the other 
hand, the cusps of the strong-cusp galaxies are ex- 
pressed by the same formula, but with 7 ~ 1 and 
a ~ 2. In addition, it was found that bright galax- 
ies have weak cusps, while faint galaxies have strong 
cusps. 

Previously, bright elliptical galaxies were believed 
to have cores with constant surface brightness. High- 
resolution surface photometry showed that these cores 
are non-isothermal ( Lauer 1985| ; Kormendy 1985) 



HST observations revealed these cores are actually 
weak cusps. 

The origin of such weak cusps is not well under- 
stood. One possible way to form a density cusp is 
to introduce the central black hole. Young (1980) 
studied the response of the structure of a stellar 
system to the adiabatic growth of the central black 
hole. He showed that the power-law density cusp with 
p oc r~ 3 / 2 is formed around the black hole. The same 
result was obtained by Quinlan, Hernquist, & Sig- 
urdsson (1995). 

Another way to form cusps using a central black 
hole was proposed by Bahcall & Wolf (1976). They 
investigated the evolution of the structure of a stellar 
system around a black hole driven by two-body relax- 
ation. They found a steady state solution with a steep 
cusp. Shapiro & Marchant (1978) and Marchant & 
Shapiro (1979, 1980) performed Monte Carlo simula- 
tion and found that this steady state cusp is actually 
formed. 

Both models predict cusps steeper than the ob- 
served ones. Therefore the origin of these weak cusps 
cannot be explained by means of these mechanisms. 

Navarro, Frenk, & White (1996) found that the 
dark matter halos formed in standard CDM cos- 



mogony is well approximated by a "universal profile" 
which approaches pa r _1 near the center. They ar- 
gued that this result might explain the origin of weak 
cusps. Recently, Fukushige & Makino (1997) per- 
formed the same calculation as Navarro et al. (1996), 
but with much higher resolution in both mass and 
space. They obtained a nearly isothermal cusp. Thus 
the cusps formed by CDM clustering scenario are also 
too steep. 

So far, the only scenario that successfully repro- 
duced the weak cusps is the merging of two galaxies 
with central black holes (Makino & Ebisuzaki 1996, 



hereafter ME). They investigated the mergings of two 
galaxies with central black holes by ./V-body simula- 
tions and found that the merger remnants have very 
shallow (p oc r~ l or shallower) density cusps in the 
central regions. This result is in good agreement with 
the HST observations. 

The physical process by which the weak cusps are 
formed in a merger remnant was not clearly explained 
in ME. It is obvious that the black holes played some 
important roles in the cusp formation, since no sig- 
nificant cusps were formed by the merger remnant of 
two galaxies without black holes. The merging of two 
galaxies with black holes involves several processes, 
such as the sinking of two black holes to the center 
of the remnant due to dynamical friction from field 
stars, the formation and hardening of a black hole 
binary and its back reaction to field stars. We need 
simpler models with which we can study the effect of 
each elementary process separately. 

In the present work, we study the sinking of a black 
hole due to dynamical friction and its back reaction to 
the distribution of field stars. We performed N-hody 
simulation of the system of spherical galaxy model 
and a black hole placed at around the half mass radius 
of the galaxy. We found the shallow density cusp with 
p ~ r -°- 5 i s formed in the center of galaxies. Thus 
we can conclude that the sinking of black hole is the 
dominant mechanism of the formation of weak cusps. 

In the next section, we describe the numerical 
methods and initial models we used. We present our 
results in and a discussion and summary in 

2. Initial Conditions and Numerical Method 
2.1. Initial Conditions 

The initial conditions we used are summarized in 
Table |l|. Column 2 shows the initial galaxy models. 
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We used King models with the central potential pa- 
rameter W c = 9, 6 and 3 as initial galaxy models. 
The number of particles in the initial galaxies (col- 
umn 3) was 32767 or 131071, and all particles have 
equal mass. We adopted the standard unit (Hcggic 
& |vlathicu 1986 ) as the system of units, in which 
G = M ga ] = 1 and E ga i — —1/4, where G is the grav- 
itational constant, M ga \ is the total mass of the initial 
galaxy and E ga \ is the total energy of the galaxy. In 
these units, the half-mass radii of the initial galax- 
ies are 0.98, 0.80 and 0.84 for the King models with 
Wc = 9, 6 and 3, respectively. The half-mass crossing 
time is 2v2 for all models. Column 4 shows the black 
hole masses. We used different black hole masses (2, 
4, 8 % of the mass of the whole system). Column 
5 shows the initial position and velocity of the black 
hole. We tried three different set of initial conditions 
of the black hole. In "off-center" runs (runs A, F, G, 
H, I and J), we initially placed a black hole in the po- 
sition of r = 1 from the center of mass of the galaxy 
(Figure [l]). In the "spiral-in" runs (run B, C and 
D), we placed a black hole in the same position as in 
the off-center runs with some initial tangential veloc- 
ity vt- We varied vt as vk, 0.5^k and 0.25t>K, where 
(circular velocity at the initial 
which is the "on-ccnter" 



r=l 

In the third set, 



v K = y/GM(r)/r 
position) . 

run (run E), however, we put the black hole exactly 
at the center of the galaxy, to compare the result with 
those of other runs. For both off-center and on-center 
runs, the initial velocity of the black hole was zero. 

2.2. Numerical Method 



We used the NBDDY1 calculation code (Aarseth 



1985|) modified for GRAPE-4, a special-p urpose corn- 



inter for astronomical ./V-body problem (Taijietal 
1996 ). In this code, we used the Hermite scheme 



( Makino fc Aarseth 1992 ) to carry out time integra- 
tion of orbits. The equation of motion for star i is 
given by 



d 2 r t 
dt 2 



§ (|r< 



Grrij (r j — r.j 



ri |2 + e 2)3/2 



(1) 



where r and m are the position and mass of the star, 
respectively, and e is the softening parameter which 
is added to avoid divergence of the gravitational force 
between close stars. The softening length is 1/1024 
for the force between field star particles and is 10 -8 
for the force between a black hole particle and a field 
star. These softening parameters are small enough 
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Fig. 1. — The snapshot of the whole system at t = 
for run A. The filled circle on (x, y) — (0, 1) shows the 
initial position of a massive BH. 

to achieve necessary spatial resolution, since we are 
interested in the central structure of the scale of ~ 
10 -2 or larger. The maximum relative energy error 
at the end of the run was 1.079 x 10~ 4 (for run A) and 
typical error was 10 -5 along the calculations. The 
black hole is represented by a particle that has much 
larger mass than that of a field particle. The mass 
accretion to the black hole through tidal disruption 
of field stars or relativistic effects were not taken into 
account. 

3. Results 

3.1. Motion of the Black Hole 

Figure || shows the trajectory of the black hole in 
the center-of-mass reference frame for the off-center 
run (run A) and the spiral- in run (run B). The time 
ranges from t = to t = 10 for run A and from 
t = to t — 15 for run B. In this coordinate 
system, the massive black hole is initially put on 
(x,y,z) = (0,1,0). In the off-center run (Figure ||a 
and c), the black hole oscillates around the center of 
the system and its amplitude is damped quickly due 
to the dynamical friction of the field stars. In the 
spiral-in run (Figure ^|b), the black hole sinks into 
the center along the spiral orbit. Even after the os- 
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Table 1 
Initial Conditions 



Run 


Galaxy Model 


N 


M BH 


BH Initial Condition 


(1) 


(2) 


(3) 


(4) 


(5) 


A 


King (Wc = 9) 


32767 


1/24 


r = 1 , Vt = 


B 


King (W c = 9) 


32767 


1/24 


r = 1, vt = v K 


C 


King (W c = 9) 


32767 


1/24 


r = 1, v t = 0.5u K 


D 


King (W c = 9) 


32767 


1/24 


r = 1, vt = 0.25«k 


E 


King (W c = 9) 


32767 


1/24 


r = 


F 


King (W c = 9) 


32767 


1/49 


r = 1, v t = 


G 


King (W c = 9) 


32767 


2/23 


r = 1, v t = 


H 


King (W c = 6) 


32767 


1/24 


r = 1, i»t = 


I 


King (W c = 3) 


32767 


1/24 


r = 1, u t = 


J 


King (W c = 9) 


131071 


1/24 


r = 1, vt = 



dilation or spiral motion is damped (t J> 5 for run A 
and t > 10 for run B), the black hole moves randomly 
around the bottom of the potential well, because the 
random velocity of the black hole reaches equiparti- 
tion with that of field stars. The typical length scale 



of this random motion is, from Figure 2|c, around 0.01 



In Figure [|, we can see that the central region of 
the initial King profile was destroyed by the black hole 
and a very weak cusp of p oc 7-~ - 5 is formed. This 
result is similar to that of M E and in good agreement 
with the HST observations ( Lauer et al. 1995 ;_ Byun 



et al. 1996 



Thus, we can analyze the system, assuming that the 
black hole is at the center, up to the radius of 0.01. 
The all following results in this paper are described 
on this assumption. The structure smaller than this 
radius is probably meaningless. 

3.2. Density Profiles 

Figures ||, ^ and || show the spatial and surface den- 
sity profiles for the off-center run (run A) , the spiral-in 
runs (runs B,C and D), and the on center run (run 



E). As mentioned in §3.1, we take the position of the 
black hole as the center of galaxy (and thus as the 
origin of the coordinate system) in the results of run 
A for t J£ 5 and runs B, C and D for t ^ 10. We 
calculated the surface density S(i?) by summing up 
the contribution of all stars outside radius R: 



£(i?) = 2 



rp 



1 

2^ 



Vr2 

E 

i(r t >R) 



zdr 



R 2 
rrij 



R? 



(2) 



where r, and rrii are the distance from the center of 
the galaxy and the mass of star i, respectively. 



Faber et al. 1996 



pebhardt et al. 1996| 
Kormcndy ct al. 1996| ). For the spiral-in runs (Figure 
, we obtained the density profiles similar to the off- 
center run. The central region has p oc r~ a5 cusp 
for all cases. On the other hand, when the black hole 
is placed initially at the center, the density structure 
becomes quite different (Figure |a and b). The cusp 
in Figure |^a is approximately p oc 



-1.5 



3.3. Velocity Structure 

In Figures |[ [?|, |s| and |9], we show the velocity pro- 
file, anisotropy and kurtosis. 

The projected (line-of-sight) velocity dispersion is 
calculated using the following formula ( Mcrrificld & 



Kent 1990) 



CTlos(-R) = 



\Z(R) 



dr rp 
,; Vr 2 - R 2 
R 2 \ . 2N IR 2 



1/2 



E 



27rS(i?) — r 

v ; i(r t >R) '" 



R 2 
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Fig. 2. — (a) Left: Trajectory of a BH on x — y plane 
for t = — 10 for off-center BH (run A); Right: Time 
evolution of the y-component of the BH's position, 
(b) Same as (a) but for spiral- in BH (run B). (c) 10 
times magnified view of the central region of (a) . 
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Fig. 3. — Results for off-center BH (run A): (a) den- 
sity profile; (b) surface density profile. 
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Fig. 4. — (a),(b) Same as Figure || but for spiral-in 
BH (run B, C and D). 
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Fig. 5. — (a),(b) Same as Figure |§| but for on-center 
BH (run E). 
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rf 



1R 2 



1/2 



,(3) 



where 



and 



are, respectively, the radial and 



tangential mean square velocities. 

In the central region (r < 0.1), the velocity disper- 
sion is a ~ r" 1 / 2 . This is simply because the grav- 
itational potential is dominated by the central black 
hole. In Figure 0, we can see that the projected one- 
dimensional velocity dispersion also shows a marked 
increase at the center. The velocity dispersion profiles 
are rather similar for off-center run and on-center run, 
since the central potential is dominated by the black 
hole in both runs. However, if we compare Figures 

and [t|c closely, it is clear that the central increase 
of (Ti os is less pronounced in the case of the off-center 
run (Figure ^a). That is because the density cusp is 



we akei. The restrits of the spiial-iii nuis (Figiui 



an d |7jb) are ahiiusL the same as LhaL of Lhe off-LeiiLei 

run, like the result for density profile. 

In Figure g, the velocity anisotropy parameter 

0=l-{v*)/2{v 2 r ). (4) 

is plotted against the distance from the center of the 
galaxy. From this equation, (3 = 0, — oo, +1 for an 
isotropic velocity distribution, purely circular orbit, 
and purely radial orbit, respectively. 

Here, the difference between off-center and on- 
center runs is quite clear. 

In run E, the change in the anisotropy is mono- 
tonic, and inner region is more radial. This can be 
interpret as the instantaneous response of stars to 
the injection of the central black hole. If the central 
potential suddenly deepens, the orbit becomes more 
radial. Note that this is quite different from the adi- 
abatic response of orbit studied by Young (1980) and 
Quinlan et al. (1995). In the case of the adiabatic 
response, orbits become more circular. 

In run A, however, the behavior of the anisotropy 
is more complex. In the outermost region, orbits are 
predominantly radial because the stars kicked out by 
the black hole populate the outer halo. The effect of 
the black hole explains the general increase in (3 in 
the outer region. 

However, in the innermost region, [3 decreased. 
There are at least two possible mechanisms for this 
decrease. First, the black hole should have ejected 
stars with radial orbit more efficiently than those with 
circular orbit. Second, the increase in mass at the cen- 
ter of the galaxy is not instantaneous as in the case 



of run E. Therefore, the response of orbits would be 
somewhere in between increase in (3 for run E and de- 
crease in j3 for the adiabatic response. If the response 
is closer to adiabatic, it is natural that (3 decreases. 

In the spiral-in runs (runs B, C and D), (3 shows 
the same tendency as that in the off-center run: (3 
increased in outer region and decreased in inner re- 
gion. The stars in central region of the spiral-in run 
are somewhat more circular than those in the same 
region of run A. Since the black hole deposited the 
angular momentum to the field stars, this result is 
quite natural. 

Figure ^ shows the distribution of kurtosis, i.e. the 
fourth moment of the line-of-sight velocity dispersion. 



We plot the dimensionlcss kurtosis k 



■ los 



>/*■ 



4 

los 



against the radius for runs A, B, C, D and E. Here 



the k urtosis is calculated as follows (Mcrrificld & Kent 



1990| ) 



'los 



)(R) 



dr rp 



r Vr 21 

2X 2 



R 2 



1 - 
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3R 4 



+3—(r 2 -R 2 )(v 2 v 2 ) + -- r 



8 r 4 
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y ' i(n>R) ' ' 



R 



„ R 2 / 9 
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R 2 )v 2 r 



R 2 



3i? 4 



.(5) 



In these figures, n = 3 corresponds to a Maxwellian 
velocity distribution. Larger values of K (k > 3) imply 
high- velocity tails, and values of k < 3 imply "boxier" 
distributions. 

In all cases, the kurtosis is always less than 3 and 
decreases toward the center, except in the very cen- 
tral region and in the outer region. This result is 
consistent with the result for an adiabatic growth BH 
model by Quinlan et al. (1995). 

3.4. Black Hole Mass & Core Radius 

Figure |l^ shows the time-averaged density profiles 
for runs A, F and G (M BH = 0.04, 0.02, 0.08M tot ) for 
t = 6 — 10. We can see that more massive black holes 
make larger cusp regions. We can also see that the 
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Fig. 6. — Velocity dispersion profiles, (a) for off- 
center BH (run A); (b) for spiral- in BH (runs B, C 
and D); (c) for on-center BH (run E). 
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Fig. 7. — Line-of-sight velocity dispersion profiles, (a) 
for off-center BH (run A); (b) for spiral- in BH (runs 
B, C and D); (c) for on-center BH (run E). 



A 

Ih 

V 



V 



<sa. 



v 

CV2 



V 



<sa. 



A 

1- 

? 

v 



V 



<sa. 



1 1 1 1 III 1 1 1 1 III 1 


Mil 




4 


1 


1 Mill 1 1 1 1 1 II 1 


Mill 


-(a) - t= o 








= (a) 


t= 




(King, W c =9) 










(King W„=9) 




t= 6 






3.5 




t= 6 








7) 

o 

b 














A 














o 


3 












? 

V 
















2.5 








— BH: r=l, v t =0 








- BH: 


r=l, v t =0 




i i i 







o 


i 


i i 






.01 .03 .1 .3 13 10 

r 



1 I III 



I I 1 1 1 



17 



-(b) t= (King - %=9) 



-G - 4=10-14 

■-»- 4=10-14 

(v t =0.5v K ) 
-*- 4=10-14 



(v t =0.25v K ) r d ; A I _ 

X A f* **\<L _ _i 



-BH: r=l, spiral— in 



j I 







.01 .03 



.3 

r 



10 



1 Mill 


linn i Mill 


-(c) 


t= 




(King, W =9)_ 




t= 6 


— BH: r=0 




i 


i 



.01 .03 



.3 

r 



10 



Fig. 8. — Profiles of anisotropy parameter for the ve- 
locity dispersion, (a) for off-center BH (run A); (b) 
for spiral-in BH (runs B, C and D); (c) for on-center 
BH (run E). Plotted points are the averaged points 
of the data in logarithmic bins for r. We put 10 bins 
Der decade. 
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Fig. 9. — Profiles of kurtosis of the line-of-sight ve- 
locity, (a) for off-center BH (run A); (b) for spiral-in 
BH (runs B, C and D); (c) for on-center BH (run E). 
Plotted points are the averaged points of the data in 
logarithmic bins for r. We put 10 bins per decade. 



slope of the cusp is almost the same for all three runs. 

Figure |Tl| a shows the relation between the black 
hole mass and the "core radius" . We adopted the defi- 
nition of the core radius by McMillan, Hut, & Makino 
(1990): 



(6) 



where rj is the distance of star i from the central black 
hole and pi is the local density at the position of star 
i. This definition gives a good estimate of the size of 
the "weak cusp" region (Makino & Ebisuzaki 1996). 
Figure |i"l| p shows the relation between the black hole 
mass and the mass of stars within the core radius. 

Figures [ll]a and b indicate that the core size shows 
almost linear dependence on the black hole mass. Of 
course, the relation is not strictly linear simply be- 
cause the initial King model had a finite core. Even 
for zero black hole mass the core size is still finite. 
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Fig. 10. — Density profiles for BH mass Mbh = 
1/49, 1/24 and 2/23 (run F, A and G), respectively. 
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Fig. 11. — (a) Relation between the BH mass Mbh 
and core radius r c for the run F, A and G, respec- 
tively, (b) Relation between the BH mass Mbh and 
core mass M c for the run F, A and G, respectively. 



3.5. Dependence on the Galaxy Model 

Figure [l^ shows the time-averaged density profile 
for t = 20 - 24 for run H and t — 40 - 44 for run I. 
We used less concentrated King model (W c — 6, 3) as 
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initial galaxy models for these runs. It has the core 
about 10 times larger than that of the model we used 
for all other runs. The result of run H (figure |l|a) 
shows that the structure of the core remains almost 
unchanged. The result of run I (figure |l2|b) shows 
that the flat core is replaced by a slight increase in 
the central density (p oc r~ 0,3 ). However, in this run, 
the initial potential energy of the black hole is not 
much larger than the bound energy of the stars in the 
core, since the core size is as large as the initial dis- 
tance of the black hole from the center of the galaxy. 
Therefore, this "cuspy" profile is likely to be formed 
through the process similar to that for run E, when 
we placed the black hole in the core. 

3.6. Effect of Two-Body Relaxation 

Figure |l3| shows the time evolution of the den- 
sity profile for runs with different number of particles 
(runs A and J). All other parameters are the same. 
In Figure |l3|a, it is clear that the central cusp be- 
comes steeper as the system evolves, while in Figure 
[l3| b (run J) the change is negligible. This difference is 
due to the difference in the thermal relaxation time. 
The model in run A has much shorter relaxation time. 
Real galaxies have the relaxation time several orders 
of magnitude larger than that of the model in run 
J. Therefore, in real ellipticals the thermal evolution 
would be negligible and the profile would not change. 

4. Discussion and Summary 

We investigated the dynamical effect of a massive 
black hole which sinks into the galactic center to the 
structure of the galaxy by means of iV-body simula- 
tions. We obtained three main results. 

First, we found that the black hole destroys the 
central core of the initial galaxy by heating up the 
field stars. After the black hole settled in the center, 
a weak density cusp (p ~ r ~ ' 5 ) remained around the 
central black hole. This result does not depend on 
how much angular momentum the sinking black hole 
initially has. This result nicely explains the recent 
HST observations of large elliptical galaxies 
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Lauer et 



Faber 



d. 1995|; [Byun ct al. 199^;|Gebhardt et al. 1996 
t al. 1996 ; Kormcndy et al. 1996j ). On the other 



hand, when a massive black hole is placed initially at 
the galactic center, the central density increases and 
a steep cusp (p oc r~ 15 ) is formed. 

Second, the velocity anisotropy profiles are quite 
different between the above two cases. For the off- 
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Fig. 12. — (a) Density profile for King model with 
W c = 6 (run H). (b) Same as (a) but for King model 
with W c = 3 (run I). 
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center black hole model, the orbits of the stars in 
the outer region tend to be radial, while in the inner 
region, the orbits tend to be circular. If the black 
hole spirals into the center, the field stars get angular 
momentum from the black hole and the orbits become 
more circular in the innermost region. For the on- 
center black hole model, however, the velocity is more 
radial near the center and approximately isotropic in 
the outer region. 

Third, the size of the weak cusp region shows a 
clear correlation with the black hole mass. The heav- 
ier black holes form larger cores. 

From the present results, we can conclude that the 
weak density cusp can be formed if a massive black 
hole dynamically sinks to the center of a galaxy. This 
suggests that the prime cause of the formation of the 
weak cusp in simulations by ME is the back reaction 
to the sinking black holes. Moreover, the origin of the 
weak cusp in large elliptical galaxies found with recent 
HST observations can be explained by this galaxy 
merging and subsequent sinking of black holes. From 
this point of view, the dichotomy of elliptical galaxies 
("core" galaxies and "power-law" galaxies (Lauer et 
al. 1995|; |gyun ct al. 1996fc |Gcbhardt et al. 1996| ; |Fabei 



ct al. 1996 ; Kormcndy ct al. 1996 )) could be under- 
stood as follows: small and therefore "cuspy" ellipti- 
cals and normal spirals are most likely to be formed 
through dissipational collapse. During the formation 
process, they are likely to contain moderately massive 
central black holes (~ 1O 7 M ). When these galaxes 
merge, black holes spiral into the center of the merger 
and left the weak cusp. If the merger still contain gas, 
significant fraction of that mass would accrete to the 



Fig. 13. — (a) Time evolution of density profile for 
N = 32767 (run A), (b) Same as (a) but for N = 
131071 (run J). 



black hole (Taniguchi & Wada 199E). In this process, 
"cuspy" profile might be restored. However, if gas is 
already consumed by some process, the weak cusp is 
preserved. Thus, massive ellipticals are likely to have 
weak cusps, while less massive ones are likely to have 
deep cusps, even if they are mergers. 

In addition, we investigated the possibilities of 
identifying the mechanisms of the formation of the 
central cusps in elliptical galaxies by comparing our 
results with the observations. We calculated the line- 
of-sight velocity dispersion and the kurtosis for both 
models, but the difference in the profiles of them 
are quite small. There was a clear difference in the 
anisotropy. However, to determine anisotropy from 
observation is not easy. Thus, it will be very diffi- 
cult to detect the difference of the two types of the 
cusps in the spectroscopic data. The most clear differ- 
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ence is in the density profile. If the improved velocity 
dispersion data for the observed "weak-cusp" galaxies 
demonstrates the presence of massive black holes, our 
scenario would be the most plausible explanation. 
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